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The  way  In  which  the  units  of  instruction  in  an  individualized 
instruction  program  are  sequenced  and  the  routing  decisions  are  made  are 
two  components  which  are  crucial  to  the  success  of  the  program.  While 
these  two  components  have  been  the  topic  of  considerable  research,  as 
yet  no  generally  accepted  procedures  have  been  developed  for  them. 

In  this  report,  a  theory  relating  performance  on  sequentially  arranged 
units  of  instruction  is  derived,  and  a  mathematical  model  for  describing 
that  relationship  is  formulated.  Procedures  for  using  the  model  to 
evaluate  sequential  relationships  and  for  making  routing  decisions  are 
described.  A  procedure  for  estimating  the  parameters  of  the  model  are 
outlined,  and  data  supporting  the  validity  of  the  model  are  presented. 

It  is  concluded  that  the  model  and  procedures  described  appear  to  be 
useful  ones,  and  that  they  appear  to  merit  continued  research  efforts 
directed  toward  their  development.  S 
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A  Latent  Trait  Hodel  for  Use  with  Sequentially 
Arranged  Units  of  Instruction 


One  of  the  fastest  growing  areas  in  the  field  of  education  today  is  the 
area  of  individualized  instruction  -  instruction  in  which  content, 
organization,  or  pacing  is  modified  for  each  individual.  Although 
individualized  instruction  comes  in  many  forms  (e.g.,  personalized  systems  of 
instruction,  computer  assisted  instruction,  individually  prescribed 
Instruction,  and  programmed  instruction),  all  of  these  forms  share  the  same 
basic  design.  All  are  basically  sequences  of  instructional  units  through 
which  subjects  are  routed  by  means  of  a  series  of  tests. 

The  way  in  which  the  units  of  instruction  are  sequenced  and  the  routing 
decisions  made  are  two  of  the  more  crucical  components  of  any  individualized 
instruction  program.  While  they  have  been  the  topic  of  considerable  research, 
as  yet  no  generally  accepted  procedures  have  been  developed  for  these 
components.  The  purpose  of  this  paper  is  to  propose  a  new  procedure  for 
developing,  evaluating,  and  implementing  routing  procedures  for  use  with 
Individualized  instruction  programs.  Specifically,  a  model  will  be  proposed 
for  describing  the  relationship  between  performance  on  sequentially  arranged 
units  of  instruction,  and  procedures  for  using  the  model  to  evaluate 
sequential  relationships  and  for  making  routing  decisions  will  be  discussed. 
Next,  a  procedure  for  estimating  the  parameters  of  the  model  will  be 
presented.  Finally,  empirical  data  will  be  analyzed  to  demonstrate  the 
validity  of  the  model.  Before  beginning  the  discussion  of  this  procedure, 
however,  some  theory  about  the  nature  of  sequential  units  of  instruction  will 
be  presented  as  a  basis  for  the  procedure. 

Sequential  Units  of  Instruction 


Underlying  Theory 

The  basic  assumption  underlying  the  sequential  arrangement  of  units  of 
instruction  is  that  performance  on  module  2  (a  unit  of  Instruction)  requires 
the  prior  knowledge  of  the  material  contained  in  module  1  (another  unit  of 
instruction).  It  might  be  true  that  all  material  in  module  1  must  be  mastered 
before  any  (or  at  least  any  appreciable  amount)  of  the  material  in  module  2 
can  be  mastered,  or  it  may  be  the  case  that  certain  sections  of  module  1  are 
prerequisite  for  certain  sections  of  module  2.  For  this  paper,  the  former 
will  be  assumed  to  be  the  case.  It  will  also  be  assumed  that  the  tests  that 
measure  the  skills  taught  in  the  modules  measure  a  unidimensional  trait. 

To  say  that  a  sequential  relationship  exists  means  that  a  certain  level 
of  performance  is  required  on  module  1  (m^)  before  learning  on  module  2  (m2) 

can  begin.  Once  that  level  (cp  is  achieved,  learning  on  mj  can  begin. 

Improvement  above  level  on  m^  facilitates  improvement  on  n^.  Once  the 

mastery  level  on  m^  (C2)  is  achieved,  additional  learning  on  m^  does  not 

facilitate  learning  on  t^.  This  relationship  is  illustrated  in  Figure  1. 


The  relationship  represented  in  the  figure  is  called  a  module  characteristic 
curve,  or  MCC. 


In  Figure  1,  the  vertical  axis  is  the  proportion  of  examinees  passing  n^, 
and  the  horisontal  axis  is  the  examinee's  status  (level  of  achievement) 

on  S|.  As  can  be  seen,  the  relation  is  horizontal  until  the  level  of 

achievement  on  m^  designated  by  c^  is  reached.  At  that  point  a  linear 

relationship  betveen  status  on  m^  and  performance  on  is  depicted.  When  the 

m^  mastery  level,  C£,  is  reached,  a  horizontal  relation  is  again  present, 

indicating  that  further  Improvement  on  m^  does  not  aid  performance  on  o^.  Of 

course,  the  relationship  in  the  range  from  c^  to  C2  need  not  be  a  linear  one, 

and  in  reality  examinees  would  be  expected  to  fall  in  a  scatter  around  the 
curve  shown  in  Figure  1. 

Figure  1 

Theoretical  Relationship 
between  Performance  on  Two  Modules 


Status  on  m 


The  low  end  of  the  curve  shown  in  Figure  1  is  not  at  zero  on  the  vertical 
axis,  nor  is  the  top  end  at  one.  It  would  be  expected  that  some  snail  portion 
of  examinees  night  pass  ^  even  with  very  little  learning  on  a^,  This  would 
be  due  to  chance  or  other  factors,  and  would  generally  be  a 

small  proportion  of  the  total  number  of  examinees.  It  would  also  be  expected 
that  some  portion  of  examinees  who  had  mastered  ttj  would  fail  m2,  simply 

because  of  failure  to  master  the  material 

not  included  in  m^. 

An  Illustration 

In  order  to  illustrate  the  processes  described  above,  simulation  data 
were  generated  according  to  the  following  process.  Item  parameters  for  the 
three-parameter  logistic  (3PL)  model  (Blrnbaum,  1968)  were  selected  for  a 
thirty  item  module  (module  2).  These  values  are  shown  in  Table  1.  Examinee 
i>i  achievement  levels  were  randomly  selected  from  a  normal  distribution  with  a 

mean  of  zero  and  a  standard  deviation  of  1.5.  The  c^  value  was  set  equal  to  6 

(achievement  level)  ”  -1.0,  and  C2  was  set  equal  to  0“  0.5.  Mastery  of  i>2  was 

arbitrarily  defined  as  seventeen  correct  out  of  thirty  items. 

For  each  examinee,  an  nu,  achievement  level  was  selected  as  follows.  If 

the  examinee's  m^  achievement  level  (0^)  was  less  than  Cj»  the  examinee's  i>2 

achievement  level  (62)  was  randomly  selected  from  a  normal  distribution  with  a 

mean  of  -1.0,  and  a  standard  deviation  of  0.5.  If  8^  >  Cj,  but  6^  <  c2 ,  62 

was  randomly  selected  from  a  normal  distribution  with  a  mean  of  6  ^  and  a 

standard  deviation  of  0.5.  If  8^  >  c^,  was  raQdo,lly  selected  from  a  normal 

distribution  with  a  mean  of  0.5  and  a  standard  deviation  of  0.5.  Table  2 
presents  a  summary  of  the  relationship  between  0^  and  Using  an<*  the 

module  2  item  parameters,  response  data  were  generated  for  module  2  according 
to  the  3PL  model  for  1000  examinees. 
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Table  2 

Summary  of  Relationship  Between  6^  and  6^ 


01 

>  C1 

02~  N(-1.0,  0.5) 

C1  <  el 

<  C2 

02~  N(0r  0.5) 

c  <  0 
c2  '  1 

02~  N(0.5,  0.5) 

Figure  2  shows  a  plot  of  6^  by  62  for  the  1000  simulated  examinees.  As 
can  be  seen,  below  -1.0  on  the  0^  scale  there  Is  a  correlation  of  about  zero 
between  0^  and  02*  Between  0^  *  -1.0  and  0^  ■  0.5  there  is  a  positive 
correlation  between  0^  and  Above  0j  *  0.5,  there  is  again  no  correlation 

between  0^  and  @2* 

Figure  3  shows  an  empirical  MCC  for  the  generated  data.  The  empirical 
MCC  was  computed  by  grouping  examinees  Into  0. 1-lntervals  of  the  ability  scale 
on  the  basis  of  0^.  For  each  interval,  the  proportion  of  examinees  in 

that  interval  who  passed  m2  was  computed  and  plotted  against  the  Interval 

midpoint.  As  can  be  seen,  the  plotted  values  form  a  rough  approximation  to 
the  curve  shown  in  Figure  1. 

The  Procedure 

The  Model 

The  procedure  proposed  for  use  with  sequential  units  of  instruction  is 
based  on  the  notion  of  the  MCC.  An  MCC  describes  the  probability  of  passing  a 
unit  of  instruction  (module)  conditional  on  latent  ability  (achievement  level) 
on  the  prerequisite  module.  The  form  of  the  MCC  proposed  in  this  paper  is  the 
four-parameter  logistic  (4PL)  model,  which  is  given  by 

"  Cj  +  (1  -  Cj  -  6j)  [l  +  EXP  (-Daj  (0lk  -  bj))]"1  ,  (1) 

where  0lk  is  the  latent  ability  of  examinee  k  on  module  i  (the  prerequisite 
module),  Pj(0^k)  i«  the  probability  of  passing  module  j  given  ability  0ik,  aj 


Is  a  slope  parameter  associated  with  module  j ,  b,  is  a  location  parameter 

associated  with  module  j ,  Cj  is  a  lower  asymptote  parameter  for 

module  j,  e^  is  an  upper  asymptote  parameter  for  module  j,  D  -  1.7,  and  EXP(x) 

■  ex.  The  Cj  term  is  used  to  account  for  the  nonzero  probability  of  passing 

module  j  for  examinees  with  very  low  ability  on  module  i,  and  the  e^  term 

accounts  for  the  nonunity  probability  of  passing  module  j  for  examinees  of 
very  high  ability  on  module  i. 

Figure  4  shows  a  4PL  MCC.  The  a-parameter  is  related  to  the  slope  of  the 
MCC  at  the  point  of  inflection,  while  the  b-parameter  serves  to  locate  the 
point  of  inflection  on  the  ability  scale. 

Using  the  Model 


Interpreting  the  Parameters.  Using  the  4PL  model  in  conjunction  with 
sequential  units  of  instruction  Involves  estimating  and  interpreting  the 
parameters  of  the  model.  The  slope  of  the  MCC,  as  indicated  by  the  a- 
parameter,  represents  the  strength  of  the  sequential  relationship.  A  steep 
slope  indicates  that  small  increases  in  achievement  on  the  prerequisite  module 
yield  large  increases  in  performance  on  the  subsequent  module.  This  would  be 
indicative  of  a  strong  sequential  relationship.  A  relatively  flat  MCC 
indicates  that  even  large  increases  in  achievement  on  the  first  module  do  not 
yield  substantial  improvement  in  performance  on  the  second  module.  This  would 
be  indicative  of  a  weak  sequential  relationship.  Thus,  the  a-parameter  serves 
as  an  indicant  of  the  strength  of  the  sequential  relationship. 

The  b-parameter  helps  to  indicate  what  level  of  performance  is  required 
on  the  first  module  to  attain  a  given  level  of  performance  on  the  second 
module.  If  the  b-parameter  for  the  MCC  shown  in  Figure  4  were  increased,  the 
curve  would  be  shifted  to  the  right.  If  this  were  the  case,  a  greater  level 
of  ability  would  be  required  on  module  1  to  attain  the  same  level  of 
performance  on  module  2  as  was  the  case  before  the  curve  was  shifted.  Thus, 
the  b-parameter  locates  the  module  on  the  achievement  scale.  The  c-parameter 
is  a  'pseudo-guessing'  parameter.  It  represents  the  probability  of  passing 
module  2  even  when  little  or  none  of  the  material  of  module  1  has  been 
mastered.  A  large  value  for  c  indicates  that  module  2  can  be  passed  at  a 
fairly  high  rate  without  knowledge  of  the  material  in  module  1.  Thus,  the  c- 

parameter  is  an  indicant  of  the  degree  to  which  module  2  can  be  passed  without 

knowledge  of  module  1  material. 

The  e-parameter  is  a  reflection  of  the  fact  that  module  2  contains 
instruction  and  material  beyond  those  in  module  1.  Perfect  mastery  of  module 
1  does  not  guarantee  mastery  of  module  2.  That  is,  module  1  is  necessary  but 
not  sufficient  for  module  2.  The  greater  the  value  of  e,  the  greater  the 
chance  of  falling  module  2  even  if  module  1  has  been  mastered. 

Setting  a  Pass/Fall  Score.  The  goal  of  setting  a  pass/fail  cut  score  for 

module  1  is  to  minimize  the  number  of  examinees  who  cannot  pass  module  2  but 
are  allowed  to  proceed  beyond  module  1  and  to  minimize  the  number  of  examinees 


Probability  of  Passing  m2 
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Figure  3 

Empirical  MCC  for  Generated 
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Figure  4 
A  4PL  MCC 


Achievement  Estimate  on  mj 


who  could  have  passed  nodule  2  but  are  held  back.  If  these  two  types  of 
errors  are  considered  equally  serious,  then  the  most  obvious  procedure  for 
setting  a  cut  score  for  module  1  is  to  determine  the  level  of  ability  on 
module  1  for  which  the  predicted  probability  of  success  on  module  2  is  0.5. 
Setting  equation  1  equal  to  0.5  and  solving  for  6  yields 


ec-  °«) +b  ■ 


where  6  is  the  pass/fail  cut  score  for  module  1 ,  ln(x)  is  the  log  to  the  base 
c 

e  of  x,  and  the  other  terms  are  as  previously  defined. 

Once  estimates  of  the  MCC  parameters  and  examinee  ability  parameters  have 

A 

been  obtained,  6  is  calculated  from  (2).  Examinees  for  whom  6  (estimated 
c 

achievement)  >  6  considered  masters  of  module  1  and  are  routed  to  module  2. 
c 

A 

Examinees  with  6  <  8  are  considered  nonmasters  and  are  not  allowed  to  proceed 

c 


to  module  2. 


Parameter  Estimation 


The  procedure  for  estimating  the  item  parameters  of  the  4PL  model 
selected  for  this  research  is  based  on  a  maximum  likelihood  estimation 
technique.  An  iterative  procedure  based  on  the  Newton-Raphson  approach  to 
solving  simultaneous  nonlinear  equations  is  employed. 

Criterion  Function 

The  estimation  procedure  is  designed  to  maximize  the  criterion  function 
given  by 


L  -  n  p,  J  Q. 

j-1  2  2 


where  L  is  the  likelihood  of  the  string  of  observed  outcomes  (passes  and 
failures)  for  a  module,  N  is  the  number  of  examinees,  Uj  is  the  module  outcome 

(zero  for  fail,  one  for  pass)  for  examinee  j,  and  Qj  is  l-Pj.  pj  is  given  by 

(1).  In  practice,  (3)  is  maximized  by  minimizing  the  negative  of  the 
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logarithm  to  the  base  e  (natural  logarithm)  of  L.  That  lst  L*  is  minimised, 

where 


* 


L 


-  loge(L) 


(4) 


Estimation  Procedure 


The  Newton-Raphson  procedure  employed  requires  the  first  and  second 
partial  derivatives  of  (4),  taken  with  respect  to  the  item  parameters.  If  f' 
is  a  column  vector  of  first  derivatives,  and  _f“  is  the  matrix  of  second 
derivatives,  then  for  any  set  of  provisional  item  parameter  estimates,  updated 
estimates  are  obtained  using  the  following  formula: 


f1+1  -  f1  +  Kf-rV]  | 

fl  » 


(5) 


where  fj"  is  the  vector  of  item  parameter  estimates  after  iteration  1,  and 
is  the  vector  of  item  parameter  estimates  after  iteration  1+1.  The  first 
and  second  derivatives  of  (4)  are  given  in  the  Appendix.  In  a  given 
iteration,  these  derivatives  are  evaluated  using  the  estimates  from  the 
previous  iteration. 

A  problem  occurs  with  the  Newton-Raphson  procedure  when  the  matrix  of 
second  derivatives,  given  by  _f" ,  is  not  positive  definite.  The  Newton-Raphson 
procedure  guarantees  convergence  only  when  _f_"  is  always  positive  definite. 

When  a  model  such  as  the  4PL  model  is  used,  the  matrix  of  second  derivatives, 
evaluated  at  the  provisional  item  parameter  estimates,  very  often  is  not 
positive  definite.  Therefore,  it  is  necessary  to  check  _f_"  for  positive 
definiteness.  If  it  is  not  positive  definite,  it  oust  be  forced  to  be 
positive  definite.  A  number  of  procedures  for  doing  this  have  been  proposed. 

Work  is  currently  underway  on  a  program  implementing  the  above  estimation 
procedure.  At  this  point  research  is  underway  to  determine  the  optimal 
procedure  for  forcing  the  matrix  of  second  derivatives  to  be  positive 
definite.  It  is  hoped  that  a  working  version  of  the  program  will  be  available 
shortly. 

Example 

In  order  to  illustrate  the  operation  of  the  estimation  procedure  just 
described,  a  preliminary  version  of  the  4 PL,  estimation  program  was  applied  to 
the  simulation  data  generated  in  the  previous  section  of  this  paper  and  for 
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which  the  empirical  MCC  is  shown  in  Figure  3.  The  true  m^  achievement  levels 
were  used  as  input  to  the  estimation  program. 

Table  3  shows  the  item  parameter  estimates  which  resulted  from  the 
application  of  the  4PL  estimation  program  to  the  simulation  data.  Figure  5 
shows  the  empirical  MCC  shown  in  Figure  3,  with  an  overlay  of  the  theoretical 
MCC  computed  using  the  item  parameter  estimates  shown  in  Table  3.  As  can  be 
seen,  the  theoretical  curve  shown  in  Figure  3  provides  a  reasonable 
description  of  the  observed  data. 


Table  3 

Item  Parameter  Estimates 
for  Simulated  4PL  Data 


Parameter 

Estimate 

a 

1.175 

b 

-0.160 

c 

0.021 

e 

0.076 

Evidence  for 

the  Validity  of  the  Model 

Method 


For  the  purposes  of  acquiring  evidence  to  either  support  or  discredit  the 
4PL  model  and  the  MCC  concept,  real  response  data  were  collected  for  a  two- 
part  arithmetic  test.  It  was  hypothesized  that  the  two  parts  of  the  test  were 
such  that  the  skills  required  for  performance  on  the  first  part  would  be 
prerequisite  to  performance  on  the  second  part.  Using  these  two  parts  as 
modules,  empirical  MCCs  were  plotted  for  various  pass/fail  cutoffs  on  the 
second  module.  These  plots  were  then  examined  as  evidence  of  the  usefulness 
of  the  4PL  model  for  use  with  these  data.  Details  of  the  process  follow. 

Data.  The  test  used  for  these  analyses  was  the  Numerical  Skills  subtest 
of  the  Career  Placement  Program  (CPP)  test  (The  American  College  Testing 
Program,  1983).  The  first  part  of  the  test,  module  1,  is  comprised  of 
nineteen  four-choice  multiple-choice  arithmetic  computation  problems,  while 
the  second  part,  module  2,  is  comprised  of  thirteen  four-choice  multiple- 
choice  word  problems  that  require  arithmetic  computation  skills  and  problem¬ 
solving  skills.  Response  data  for  these  items  were  collected  for  3768  cases 
from  the  1983  norming  administrations  of  the  test. 
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Since  there  Is  no  already  determined  pass/fail  cutoffs  for  the  CPP 
subteats,  the  analyses  performed  in  this  stage  of  the  research  were  repeated 
for  a  number  of  different  cutoffs  for  module  2,  so  as  to  avoid  any 
capitalisation  on  chance  from  the  cutoff  selection.  Using  a  given  pass/fail 
cutoff  for  module  2,  each  examinee  was  assigned  a  score  of  0  (fail)  or  1 
(pass)  depending  on  whether  the  examinee's  raw  score  on  module  2  exceeded  the 
cutoff  for  module  2.  These  0,  1  data,  along  with  examinees'  achievement  level 
estimates  from  module  1,  were  the  input  for  these  analyses. 

Ability  Estimation.  The  achievement  level  estimates  on  module  1  for  the 
examinees  were  obtained  through  the  application  of  the  3PL  model  to  the 
examinees'  response  data  for  module  1.  The  LOGIST  estimation  program 
(Wlngersky,  Barton,  and  Lord,  1982)  was  used  to  estimate  the  parameters  of  the 
3PL  model. 

Plotting  HCCs .  The  initial  step  in  these  analyses  was  the  division  of 
the  achievement  scale  into  a  number  of  narrow  Intervals  (0.1  width). 

Examinees  were  then  sorted  into  these  intervals  on  the  basis  of  their  module  1 
achievement  level  estimates.  For  a  given  module  2  pass/fail  cutoff,  the 
proportion  of  examinees  within  each  Interval  passing  module  2  was  computed. 

For  each  module  2  cutoff,  the  proportions  passing  module  2  were  plotted 
against  the  Interval  midpoints,  thus  forming  an  empirically  derived  MCC. 
Adjacent  intervals  were  collapsed  to  assure  an  Interval  sample  size  of  at 
least  ten.  These  HCCs  were  examined  to  assess  the  reasonableness  of  the  4PL 
model  for  describing  the  form  of  the  resulting  curve. 

Results 


Figures  6  through  12  show  the  empirical  HCCs  obtained  for  the  CPP  data 
for  pass/fall  cutoffs  on  module  2  of  three  through  nine  correct  out  of  the 
thirteen  items,  respectively.  Table  4  shows  the  obtained  proportions  passing 
plotted  in  Figures  6  through  12.  Table  4  also  shows  the  numbers  of  examinees 
in  the  different  Intervals. 

As  can  be  seen  from  these  figures,  the  relationship  between  module  1 
ability  and  module  2  performance  does  appear  to  be  at  least  a  monotonlcally 
increasing  one.  Also,  for  several  of  the  plots,  there  appears  to  be  a  non- 
unity  upper  asymptote.  It  is,  however,  difficult  to  discern  a  nonzero  lower 
asymptote  in  these  plots.  Of  course,  a  lower  asymptote  of  zero  is  a  special 
case  of  the  4PL  model.  It  may  eventually  be  fruitful  to  drop  the  lower 
asymptote,  but  as  yet  there  is  little  evidence  to  support  such  a  step. 

There  are  some  interesting  trends  evident  in  Figures  6  through  12.  As 
the  pass/fall  cutoff  score  on  module  2  increases,  of  course,  fewer  examinees 
of  low  achievement  level  on  module  1  pass  module  2.  If  the  material  in  module 
2  requires  the  knowledge  of  module  1  material,  clearly  requiring  more  module  2 
material  for  passing  will  require  more  module  1  material. 

As  the  module  2  pass/ fail  cutoff  Increases,  the  upper  asymptote  of  the 
MCC  decreases  (the  e  term  increases  in  value).  This  is  an  indication  that 
complete  knowledge  of  module  1  is  not  sufficient  for  guaranteed  success  on 
module  2.  Another  way  of  saying  this  is  that  word  problems  require  more 
knowledge  than  simply  mastering  arithmetic  operations. 
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The  patterns  evident  In  these  figures  suggest  that,  if  nodule  2  were 
still  easier  to  pass  than  was  the  case  with  the  pass/fall  score  of  3,  there 
would  be  a  nonzero  lower  asymptote  to  the  MCC.  Unfortunately,  for  this 
particular  test  lower  cutoffs  yielded  an  alnost  flat  MCC  near  unity.  Almost 
all  exaninees  got  at  least  two  items  correct  on  module  2,  regardless  of  their 
module  1  ability. 


Summary  and  Conclusions 

While  this  research  project  is  still  incomplete,  it  has  yielded 
encouraging  results.  A  theory  relating  performance  on  sequentially  arranged 
units  of  instruction  was  derived,  and  a  model  for  describing  that  relationship 
was  formulated.  Procedures  for  using  the  model  to  evaluate  sequential 
relationships  and  for  making  routing  decisions  were  described.  A  procedure 
for  estimating  the  parameters  of  the  model  was  outlined,  and  data  supporting 
the  validity  of  the  model  were  presented.  All  things  considered,  the  model 
and  procedures  described  appear  to  be  useful  ones,  and  they  appear  to  meric 
continued  research  efforts  directed  toward  their  development. 
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where  Pj  ■  Pj(0^j),  Q j  -  1  -  P j ,  Xj  -  Da(8j  -  b),  and  the  remaining  terms  are 
as  defined  for  (1).  The  matrix  of  second  derivatives,  denoted  byjT',  is  given 
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The  matrix  is  symmetric.  The  individual  terms  in  the  matrix  are  given  by: 
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